A Laboratory Search for a Long-Range T-odd, P-odd Interaction from Axion-Like Particles 
using Dual Species Nuclear Magnetic Resonance with Polarized 129 Xe and 131 Xe Gas 
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Various theories beyond the Standard Model predict new particles with masses in the sub-eV range 
with very weak couplings to ordinary matter. A new P-odd and T-odd interaction between polarized 
and unpolarized nucleons proportional to s • r is one such possibility, where r is the distance between 
the nucleons and s is the spin of the polarized nucleon. Such an interaction involving a scalar coupling 
g 3 at one vertex and a pseudoscalar coupling g p at the polarized nucleon vertex can be induced by 
the exchange of spin-0 bosons. We used the NMR cell test station at Northrop Grumman Corporation 
to search for NMR frequency shifts in polarized 129 Xe and 131 Xe when a non-magnetic zirconia rod 
is moved near the NMR cell. Long (T 2 ~ 20 s) spin-relaxation times allow precision measurements 
of the NMR frequency ratios, which are insensitive to magnetic field fluctuations. Using existing 
theoretical calculations of the neutron spin contribution to the nuclear angular momentum in xenon 
nuclei, the measurements improve the the laboratory upper bound on the product g s g% by 2 orders of 
magnitude for distances near 1 mm. The sensitivity of this technique can be increased by at least two 
more orders of magnitude. 

PACS numbers: 14.20.Dh, 13.75.Cs, 14.80.Va, 24.80.+y 



The possible existence of new interactions of nature 
with ranges of macroscopic scale (millimeters to mi- 
crons) and very weak couplings to matter was sug- 
gested long ago Hi 0] and has recently begun to attract 
more scientific attention. Particles which might trans- 
mit such interactions are starting to be referred to gener- 
ically as WISPs (Weakly-Interacting sub-eV Particles) |3|] 
in recent theoretical literature. Many theories be- 
yond the Standard Model possess extended symmetries 
which, when broken at a high energy scale, produce 
pseudo-Nambu-Goldstone bosons and lead to weakly- 
coupled light particles with relatively long-range inter- 
actions such as axions, familons, and Majorons [4]. Sev- 
eral theoretical attempts to explain dark matter and dark 
energy also produce new weakly-coupled long-range 
interactions. The fact that the dark energy density of or- 
der (1 meV) 4 corresponds to a length scale of 100 fxxn 
also encourages searches for new phenomena around 
this scale j7t|j|]. Long ago Moody and Wilczek 0] con- 
sidered the form of interactions which could be induced 
by the exchange of a spin field between fermions 
with scalar or pseudoscalar couplings. They highlighted 
an interesting P and T violating scalar-pseudoscalar 
(monopole-dipole) interaction of the form 



V = h 2 g s g p - 



a ■ r 



87rm„ V ?*A r 2 



(1) 



where m„ is the mass and s = Ha/2 is the spin of 
the polarized particle, A is the interaction range, f = 
r/r is the unit vector between the particles, and g s 
and g p are the scalar and pseudoscalar coupling con- 
stants. Axions |Tj3l . axion-like particles |0] predicted 
by string theory fjj, or a very light spin-1 boson (lH) . 
can induce such an interaction and are candidates for 
cold dark matter The axion mass is constrained 

to the so-called "axion window" 0] between 1 fxeV 
to 1 meV, corresponding to a force range between 2 
cm to 20 fim. Most of the experiments which have 
been performed to search for such interactions il3Ml8h 
are sensitive to ranges A > 1 cm. Recent experimen- 
tal work to constrain the monopole-dipole interaction 
at shorter distances has employed ultracold neutron 
bound states I19l l20l] . polarized ultracold neutrons in 
material t rap s | 2lL |22|], and spin relaxation of polarized 
nuclei ESS], with the most stringent direct experimen- 
tal constraint to date on the coupling constant product 
g s g p coming from an analysis of NMR measurements 
on polarized 3 He |2o1. |2^| . Torsion balance measure- 
ments have recently set new stringent limits on possible 
monopole-dipole interactions involving polarized elec- 
trons 12811 . Although all of these laboratory results are 
less stringent than those inferred from separate limits 
on g s and g p through a combination of torsion balance 
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experiments and astrophysical constraints from stellar 
evolution and SN1987A 12911 , direct laboratory experi- 
ments are more direct and therefore are of fundamental 
interest. We used a NMR cell test station at Northrop 
Grumman Corporation (NGC), originally developed to 
evaluate dual species xenon vapor cells for a nuclear 
magnetic resonance gyroscope, to search for NMR fre- 
quency shifts in polarized 129 Xe and 131 Xe when the 
position of a non-magnetic zirconia rod is modulated 
near the NMR cell. By comparing the NMR frequen- 
cies of the two Xe isotopes, magnetic field changes can 
be distinguished from frequency shifts due to the pos- 
sible monopole-dipole coupling of the polarized Xe nu- 
clei to the nuclei in the zirconia rod. The simultaneous 
observation of precession of both isotopes allows us to 
eliminate possible systematic effects from magnetic field 
fluctuations. We use recent theoretical calculations lf30h 
of the neutron spin contribution to the nuclear angular 
momentum in 129 Xe ( J = 1 /2) and 131 Xe (J = 3/2), orig- 
inally conducted to evaluate their sensitivity for spin- 
dependent scattering of dark matter, to put a new upper 
bound on the product g s g™ for ranges at the mm scale. 

The measurements were performed using a magneti- 
cally shielded co-magnetometer cell consisting of 85 Rb, 
5 Torr 129 Xe, 45 Torr 131 Xe, and 250 Torr N 2 . Figure [T] 
shows a schematic of the apparatus. The Rb atoms 
are spin-polarized by absorption of circularly polarized 
light from a pump laser and, through spin-exchange 
collisions, in turn polarize the Xe nuclei parallel to a 
DC magnetic field [31]. The Rb atoms also serve as a 
sensitive high-bandwidth magnetometer that simulta- 
neously detects the Larmor precession of the two xenon 
isotopes; the transverse magnetic fields of the polarized 
nuclei produce an oscillating transverse spin polariza- 
tion of the Rb atoms that is detected as a rotation of the 
polarization of the linearly polarized sense laser. 

We record simultaneous free-induction-decay (FID) 
signals from both 129 Xe and 131 Xe . We initiate the FIDs 
by applying short resonant ~ ir/4 pulses to rotate the 
two Xe polarizations away from the z-direction. We iso- 
late a single species FID by multiplying the raw mag- 
netometer data, which is a superposition of the FIDs 
for both isotopes, by cos(27r/ i) where /o is chosen to 
be typically 0.5 Hz below the resonance frequency for 
the species of interest. We then take a FFT of the data, 
apply a 20 Hz low-pass filter, and inverse transform 
to get the frequency down-shifted FID for the species 
of interest. We eliminate the first 1.6 seconds and the 
last 0.2 seconds of this signal, which contain filter arti- 
facts and pulse residuals. Over the 100 sec FID acqui- 
sition time, there are Zeeman frequency drifts that may 
arise from alkali polarization drifts or actual magnetic 
field drifts. We take the 129 Xe signal as our magnetome- 
ter for isolating these drifts, fitting the 129 Xe FIDs to 
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FIG. 1: Dual species free-induction-decay spectrometer. A 2 
mm glass cell containing 129 Xe and 131 Xe polarized by spin- 
exchange optical pumping of 85 Rb is in a uniform magnetic 
field inside a magnetic shield. A moveable zirconia rod enters 
from one end of the apparatus and the pump and probe lasers 
enter from the opposite side. Free induction decay signals of 
the two Xe isotopes are sensed optically using the Rb vapor 
as a magnetometer. Sample FIDs, with fits, show 129 Xe above, 
131 Xe below. 



the form Si = Acos(2nfi(t + p(t)) + 0i)e~ t/T21 where 
p = SB/B parametrizes the fractional magnetic field 
drift with time. We use the deduced p(t) as input for 
the fits to the 131 Xe FIDs to deduce the 131 Xe frequency 
/2 as described in more detail below. Sample FIDs are 
shown in Fig.[TJ 

We performed the g s g™ search by comparing the ra- 
tios /1//2 when a mass was close to or far away from 
the cell. To this end we attached a mounting structure to 
the magnetic shielding ( DC shielding factor ps 3 x 10 4 ) to 
allow the repeatable positioning of a zirconia rod close 
(750 ± 60 microns from the inner surface) far (1.0 ± 0.1 
cm) from the inner surface of the cell. The vapor cell 
was a 1.9 mm internal dimension cube with a 0.5 mm 
diameter fill stem pinch-off on one side and nominally 
0.55 mm thick walls. The rod position was alternated 
for successive runs, with the time from the start of one 
data collection run to the start of the next being approx- 
imately two minutes (4 — 5 times T\ for each xenon iso- 
tope). Either the polarization of the pumping light or 
the direction of the magnetic field was reversed every 
10th run, and the magnitude of the magnetic field was 
changed once, from 0.13 G to 0.064 G. 

The J = 1/2 129 Xe isotope couples only to vector 
fields such as the external magnetic field, whereas the 



J = 3/2 131 Xe isotope can also experience quadrupole 
interactions with the electric field gradients on the cell 
walls. The nuclear spin Hamiltonian for isotope j is 

Hj = h l3 B K z + hXjK z + S jt2 H Q (2) 

where j = 1 for 129 Xe and j = 2 for 131 Xe. Thus 
there are 4 NMR frequencies deduced from the FIDs: 
h = \lx\B! + z ■ BXi, hi = \l2\B2 - z ■ BX 2 + q. t 
where the gyromagnetic ratios of species j are 7j, the 
monopole-dipole couplings are X,, and the quadrupole 
shifts of the 3 NMR lines of 131 Xe are qi (i = -1, 0, 1). 
The presence of the three frequencies is evident in the 
unusual shape of the 131 Xe FID in Fig. Q] We allow for 
the possibility that the two Xe isotopes may experience 
slightly different magnetic fields Bj due to slightly dif- 
ferent interactions with the polarized Rb atoms. In pre- 
vious experimental work using J = 3/2 nuclear po- 
larized species to search for various exotic physical ef- 
fects IB2I434I1 , the quadrupole interactions were typically 
treated in a spherical approximation. Our cell is nomi- 
nally cubic but the 6 internal surfaces are not identical, 
and the glass stem is located on one of the cell faces 
whose normal is along the y-direction. The averaged 
quadrupole interaction on a surface with normal h is 
V n = Qn (K • nn • K - K ■ K/3) /2 We therefore 

consider the quadrupole interaction as a sum of contri- 
butions from each of the 3 pairs of orthogonal faces of 
the cell. 2nd order perturbation theory shows that the 
three resonance frequencies in the free-induction decay 
are 

hi = h2B 2 \-z-BX 2 + 2i\Q ll \ + ^^ (3) 

where Q\\ =Q Z - (Q x + Q y )/2 and Q ± =Q X - Q y . The 
first and second order quadrupole shifts are typically 20 
mHz and 500 /iHz respectively. By taking the average 

h = hl Y 2 ' 1 = -z-BX 2 (4) 

the quadrupole interaction is cancelled to second order. 

The magnetic field experienced by the two nuclei is 
not only the applied field, but contains an additional 
contribution B\jP from the alkali atoms of polarization 
P[36]. This alkali field is about 0.1 mG in this experi- 
ment. To the extent that the alkali field is independent of 
the Xe isotope, the noble gas frequency shift will be pro- 
portional to the nuclear gyromagnetic ratio. As shown 
in Fig. [2j we observe a small shift in the frequency ra- 
tio that is correlated with the product of the sign of 
the light polarization and the magnetic field direction. 
This signals the presence of a slight isotopic difference 
in the alkali fields experienced by the two nuclei, to our 
knowledge the first observation of this phenomenon in 
spin exchange optical pumping. We account for this 
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FIG. 2: A sampling of frequency ratio measurements, showing 
a correlation with the product of the magnetic field direction 
and the polarization direction, as expected for a differential 
alkali field shift. 



by parameterizing it as Bai = Ba(1 + SBa/2), Ba 2 = 
Ba(1 — SBa/2). To first order in the small shifts v ro t, 
SB a = ABa/Ba, and X, the frequency ratio is 

— = r (1 + SBa—s— )+z-B (5) 

/2 V B J /2 

where the ratio of gyromagnetic ratios is r = |ti | / |T2 1 - 
A monopole-dipole coupling causes a rod in/rod out 
change in /1 / f 2 that reverses with magnetic field direc- 
tion and is independent of polarization. The isotope 
effect, on the other hand, correlates with the product of 
the magnetic field direction and the polarization direc- 
tion. In Fig. 12 the magnetic field magnitude was halved 
at run 600, causing a doubling in the modulation ampli- 
tude that is reproduced in the data, consistent with Eq.[5j 
The fractional isotope shift in the alkali field, deduced 
only from runs with the rod out, is SB a = 0.0017(1). 

The whole sequence of frequency shift measurements, 
corrected for the differential alkali field shift, are shown 
in Fig. [3] The data were taken over 4 days, with vary- 
ing degrees of magnetic noise and optimization of mag- 
netometer signals accounting for variations in the noise 
levels seen. About 20 runs were discarded due to un- 
usual magnetic field variations that our fitting was not 
able to account for. Another set of runs during which 
the alkali field shift was rapidly varying, probably due 
to laser instabilities, were also discarded. We fit the cor- 
rected frequency ratio data to the function r + z ■ BXj j 2 
with the ratio of gyromagnetic ratios r and the weighted 
difference of monopole-dipole shifts X = X\ — rX 2 . The 
result is r = 3.3734081(5), and X = -8(23) /iHz. The 
value of r is consistent with Brinkmann 13711 , who found 
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FIG. 3: Frequency ratio measurements obtained with a variety 
of magnetic field and polarization directions, and two mag- 
netic field strengths, with and without a nearby zirconia rod. 
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FIG. 4: Constraints on the coupling constants g s g% of the spin- 
dependent force as a function of the force range A. The dashed 
line is the result from ffll . the dash-dotted line is from l26ll , 
the red line is from |2^1 , and the solid line is the present work. 



3.37340(4) and observed a pressure dependence. We 
measured only one cell and have no information about 
pressure effects, so our value of r is valid for the param- 
eters of our cell. Our number may also be affected by an 
imperfect model for the quadrupole interaction. 

To use our result for X to set a limit on the monopole- 
dipole interaction, which couples to the nucleon spins 
in the nucleus, we need information on the neutron and 
proton spin contributions to the total nuclear angular 
momentum of 129 Xe (J = 1/2) and 131 Xe (J = 3/2). 
Theoretical calculations have been done by several nu- 
clear theorists over the last decade to calculate the sen- 
sitivity of experimental searches for dark matter using 
liquid xenon to dark matter /xenon spin flip scattering. 
From these calculations |30j I38H40I1 we use the most re- 
cent effective field theory-based results of Menendez et 
al. |30|l which reproduce well the low-lying energy lev- 
els of 129 Xe and 131 Xe. According to this calculation 
(S™ 9 ) = 0.33, (S 1 , 29 ) = 0.01, (S™ 1 ) = -0.27, and 
(S 1 , 31 ) = -0.01. Thus 66% of the J = 1/2 of the 129 Xe 
nucleus comes from the neutron spin and 18% of the 
J = 3/2 polarization of the 131 Xe nucleus comes from 
the neutron spin: almost all the rest is neutron orbital 
angular momentum. Furthermore, in both 129 Xe and 
131 Xe the proton spins make almost no contribution to 
J. Therefore to an excellent approximation the results 
can be interpreted in terms of neutron couplings. 

By integrating Equation [I] over the polarized xenon 
and unpolarized zirconia mass distributions (the nu- 
cleon density of the zirconia is 3.6 x 10 24 cm" 3 ), we ob- 



tain a limit on the monopole-dipole coupling strength 
as a function of the range A which excludes an area in 
the (g s g p , A) plane above the line defined in Figure |U 
The limits shown in Figure |4] involve the pseudoscalar 
coupling g'p of the neutron and the scalar coupling g s of 
an ensemble of unpolarized matter made of a roughly 
equal proportion of neutrons and protons. 

We have set a new upper bound on the product of 
the scalar and pseudoscalar coupling constants g s g™ be- 
tween nucleons and a polarized neutron for ranges near 
mm distances. Our bounds g s g p < 10 -19 for A = 0.1 
mm to g s g p < 10~ 23 for A = 1 cm are the most stringent 
laboratory limit over these distances, which correspond 
to exchange boson masses from 10 -3 to 10 -5 eV/ c 2 . 

There is great room for improvement of the mea- 
surement sensitivity using this technique. Direct sta- 
bilization of the 129 Xe frequency using magnetic field 
feedback, plus stabilization of the alkali magnetization, 
would greatly reduce the leading uncertainties and im- 
prove the sensitivity to much less than 1 /^Hz. 
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